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A harmonically trapped ultracold ID spin-1 Bose gas with strongly repulsive or attractive ID even- 
wave interactions induced by a 3D Feshbach resonance is studied. The exact ground state, a hybrid 
of Tonks-Girardeau (TG) and ideal Fermi gases, is constructed in the TG limit of infinite even-wave 
repulsion by a spinor Fermi-Bose mapping to a spinless ideal Fermi gas. It is then shown that in the 
limit of infinite even-wave attraction this same state remains an exact many-body eigenstate, now 
highly excited relative to the collapsed generalized McGuire cluster ground state, showing that the 
hybrid TG state is completely stable against collapse to this cluster ground state under a sudden 
switch from infinite repulsion to infinite attraction. It is shown to be the TG limit of a hybrid super 
Tonks-Girardeau (STG) state which is metastable under a sudden switch from finite but very strong 
repulsion to finite but very strong attraction. It should be possible to create it experimentally by a 
sudden switch from strongly repulsive to strongly attractive interaction, as in the recent Innsbruck 
experiment on a spin-polarized bosonic STG gas. In the case of strong attraction there should also 
exist another STG state of much lower energy, consisting of strongly bound dimers, a bosonic analog 
of a recently predicted STG state which is an ultracold gas of strongly bound bosonic dimers of 
fermionic atoms, but it is shown that this STG state cannot be created by such a switch from strong 
repulsion to strong attraction. 

PACS numbers: 03.75.-b,67.85.-d 



If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excitation 
quantum is larger than available longitudinal zero point 
and thermal energies, the effective dynamics becomes 
one-dimensional (ID) [l], 3D Feshbach resonances 
Q allow tuning to the neighborhood of ID confinement- 
induced resonances 0, 0] where the ID interaction is very 
strong, leading to strong short-range correlations, break- 
down of effective-field theories, and emergence of highly- 
correlated TV-body ground states. In the case of spin- 
less or spin-polarized bosons with zero-range Lieb-Liniger 
(LL) [H| delta function repulsion gsS^Xj — xg) with cou- 
pling constant qb — > +oo, the Tonks-Girardeau (TG) 
gas, the exact iV-body ground state was determined in 
1960_by a Fermi-Bose (FB) mapping to an ideal Fermi 
gas [6j , leading to "fermionization" of many properties of 
this Bose system, as recently confirmed experimentally 
0,1- 



It was predicted several years ago 0-11 1 that if gs is 
large and negative, there exists a highly excited gaseous 
state known as the super Tonks-Girardeau (STG) gas 
which is metastable against collapse to the McGuire clus- 
ter ground state [ijj], and this state was recently cre- 
ated by the Innsbruck group [13J by suddenly switching 
the interaction from strongly repulsive to strongly at- 
tractive by passing through the ID confinement-induced 
resonance induced by a 3D s-wave Feshbach resonance. 
This initially surprising metastability was recently shown 
14| to result from the fact that after a sudden switch 
from infinite repulsion (gs = +oo) to infinite attrac- 
tion (gs = -co), the gs = +oo TG ground state re- 
mains an exact energy eigenstate, now highly excited 



relative to the collapsed McGuire cluster ground state, 
and is therefore stationary and hence completely stable 
against collapse. Then by continuity, the system should 
be metastable under a sudden switch from strong repul- 
sion to strong attraction as in the Innsbruck experiment 

A generalization of the FB mapping to a ID spinor 
Fermi gas is relevant to recent predictions of sup er Tonks- 
Girardeau (STG) states in such a system 13, LLfl] . The 
model consists of an ultracold ID gas of fermionic atoms 
with very strong attractive interactions in an effectively 
ID trap, in two different hyperfmc states. Similar STG 
states should also occur in a ID spinor Bose gas con- 
sisting of bosonic atoms in three hyperfine states conve- 
niently labeled by spin z-component labels a = — 1, 0, 1 of 
spin-1 atoms, but have not been studied so far, although 
the case of strong repulsion (TG regime) has been studied 
in detail by a spinor generalization of the FB mapping 
171 ]. The Hamiltonian used there included the kinetic 
energy, a spin-independent harmonic trapping potential, 
a zero-range LL interaction [Bj], and a spin-spin interac- 
tion, but simplifies in the TG regime where the spin-spin 
interaction is negligible. Then the Hamiltonian is of the 
same form as that in 15, [l6j], consisting of the kinetic 
energy, a spin-independent harmonic trap potential, and 
a spin-independent LL interaction: 



h 2 d 2 

2m dx 2 
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i<j<e<N 



5(xj-xt) 
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The case relevant to generation of a STG state is that 
of strong s-wave scattering due to a 3D s-wave Feshbach 
resonance, which leads in ID to a confinement-induced 



2 



even- wave resonance [l|, 0| and a LL delta function inter- 
action It acts between atoms j and £ only if the rel- 
ative wave function is symmetric under spatial exchange 
(xj <-> a;^), in which case bosonic symmetry under com- 
bined space-spin exchange (xj,aj) -H> {xt,<jg) demands 
that the relative wave function be symmetric under spin 
exchange Uj 4-> ag. On the other hand, if the relative 
wave function is antisymmetric under spatial exchange, 
then the interaction term in Eq. (TTJ) is cancelled by the 
node at Xj = xg. 

If gB is large and negative, two quite different STG 
states can occur, as in [16j. One is a hybrid between an 
STG gas with strong attractions between atoms in spin- 
symmetric pair states and an ideal Fermi gas with no 
interactions, the STG component being an exact analog 
of the STG state in an ultracold ID Bose gas predicted in 
[§- 11 1 and recently created by the Innsbruck group [l3[ . 
It should be possible to create it experimentally by a sud- 
den switch of the interaction from strongly repulsive to 
strongly attractive, as in [lj|. The other STG state is a 
trapped analog of a recently predicted STG state in the 
spinor Fermi gas [IH , an ultracold gas of strongly bound 
boson dimers which behave as bosons with a strongly 
attractive boson-boson interaction, leading to STG be- 
havior. 

Consider first the case N = 2. The analysis closely par- 
allels that for the spinless Bose gas [14[ and the spinor 
Fermi gas 16]. The LL interaction is gsS(xi — X2), the 
harmonic trap potential is muj wave func- 

tion for the center of mass (cm.) coordinate X = (xi + 
x 2 )/2 is V^.m. = exp[-X 2 /xl sc ] where x osc = yjh/mui, 
and its energy is E c . m _ = hu/2 assuming that the cm. 
mode is unexcitcd. Bosonic symmetry requires that the 
relative wave function ^rel(^i) a i] x 2> 02) be symmetric 
under conbined space-spin exchange. Since the Hamilto- 
nian is spin-independent, ip rc i factorizes into the product 
of a spatial wave function 4>(x), with x — x\ — X2, and 
a spin-pair wave function x( £r i) cr 2) I where either both 
<f> and x are symmetric or else both are antisymmet- 
ric. If they are antisymmetric then <p vanishes at con- 
tact x = 0, killing the LL interaction gB&{x\ — X2), so 
that <j) reduces to a trapped harmonic oscillator eigen- 
state odd in x, the product of the ground state Gaus- 
sian and a Hermite polynomial of odd order. In the 
symmetric case <fi is identical with the solution for spin- 
less bosons [3 HH, HH, El- These solutions are an- 
alytic continuations of the Hermitc-Gaussians to non- 
integral quantum number v, and are parabolic cylin- 
der function D v {q) where x — qx osc and the allowed 
values of v are solutions of a transcendental equation 
r(| — hv)/T(— ~v) — —A where A is a dimensionless cou- 
pling constant A = gs / (2 3 / 2 hwx asc ) [18|. This gives the 
solution for q > 0, while for q < it is D u (\q\), since 
4>{x) is even. For A^ — > +00 the solution reduces to the 
first excited harmonic oscillator state 4>{x) = |g|e _9 I 2 



[20|| , the usual TG gas ground state with a cusp at q = 
due to the point hard core interaction. For A — > — 00 
this wave function is still an exact energy eigenstatc, 
but it is highly excited, the ground state being a col- 
lapsed state which is an analog, for the trapped system, 
of McGuire's cluster state [12J. It is an even solution 
also expressible in terms of a Z?„, but one whose energy 
approaches —00 as gs — > —00 (am — > 0+); see Fig. 
5 of 18(. For aiD — > 0+ it is well approximated by 
ipso ~ exp(-|xi - x 2 \/a 1 n) exp[—(xf + x\)/2x 2 sc ] where 
Aid > is the ID scattering length, am vanishes as 
|A| — > 00 and increases with decreasing |A| [l|, [H, H], 
and at the same time the node moves from the origin 
to a position x no ^ c which is close to <Zid when |A| ^> 1 
but slightly smaller, the explicit expression being given 
in Q. 

An explicit formula for an exact N = 2 ground state 
in the TG limit Xb — > +00 can be given in terms 
of a mapping function a{xi,<j\;x 2 ,a 2 ) which is space- 
spin antisymmetric and everywhere ±1, which generates 
the strongly-interacting ground state ipso from the spin- 
less ideal Fermi gas ground state '(/'ideal, = H\{qi — 
q 2 )e~( Xl+X2%, /( 2x °'<:) where H\ is the Hermite polynomial 
of order one, the lowest antisymmetric state. The desired 
representation is ipso = a(xi, o\\ x 2 , o^V'idcaio^i, x 2 ) 
where 

1 

a(xi,(Ti ; x 2 ,a 2 ) = ^ (5<7 1 , s i5 (r2 , 3 sgn(xi - x 2 ) 

s=-l 

~\~ ^ ' (<^CTi ,Si ^0" 2 ,S2 ^<Tl ,S2 ^(T2.Sl ) ■ (2) 

-1<S1<S 2 <1 

The signum function sgn(xi — X2) introduces the re- 
quired cusps at contact x\ — x 2 due to the zero-diameter 
hard-core interaction for space-even, spin-even scatter- 
ing, which are absent in the space-odd, spin-odd part of 
the wave function; note that sgn(xi — a^^idcai.o^ii £2) 
is space-even and V'ideai.oC^ij^) is space-odd. ipso is a 
hybrid TG-idcal Fermi gas state. Since the Hamiltoman 
is spin-independent, the ground state is 3 2 = 9-fold de- 
generate, so this ipBo is only one member of the degener- 
ate ground manifold. One can obtain others by changing 
the signs of various terms in |2[). The arbitrary over- 
all sign can be determined by the s = — 1 term in the 
first summation, and changing the sign of either of the 
two others gives two additional states, for a total of 3. 
Then changing the sign of any of the three terms in the 
second summation multiplies this by 3, giving an overall 
total of 9. These are linearly independent although not 
orthonormal, and span the 9-dimensional ground mani- 
fold. In the absence of spin-spin interaction or an exter- 
nal magnetic field, these states are experimentally indis- 
tinguishable; they are chosen here for easy generation by 
mapping from the spinless ideal Fermi gas, rather than 
by requiring them to be simultaneous eigenstates of S* 2 
and S z . 
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Generalization to N > 2 is straightforward. Denote 
the iV-atom wave functions by ipB(x%, ar, ■ ■ ■ ; x^j cjv)- 
Consider first the case A — > +00. This can be gener- 
ated from an ./V-particle ID spinless ideal Fermi gas state 
V'idcai by a space-spin mapping function M which is the 
product of two-particle mapping functions @ over all 
pairs: 



ip B = M(x 1 ,<r 1 ; ■ ■ ■ ;x N ,a N )t/j:_ 



ideal 



(3) 



where the spin-dependent Fermi-Bose mapping function 
M is 

M{xx,ai\ ■ ■ ■ ;x N ,a N ) = ]J[ a(xj , aj; x e , a e ) (4) 

i<j<e<N 

and 7/>idoai is an energy eigenstate of the trapped ID 
ideal gas of "spinless" fermions, a Slater determinant of 
TV different harmonic oscillator orbitals. Starting from 
this ideal Fermi gas, this mapping generates the required 
TG-gas cusps in space-symmetric, spin-symmetric scat- 
tering channels of the strongly interacting spinor Bose 
gas, but no interaction in space-antisymmetric, spin- 
antisymmetric channels. M is space-spin antisymmetric 
2l| and everywhere ±1, hence self-inverse. The Hamil- 



tonian commutes with M when all \xji\ > where 
Xj£ = Xj — X£, and since the interaction term in ([1]) 
vanishes there ifjB is an energy eigenstate with the same 
eigenvalue as V'ideai- Here we are particularly interested 
in the ground state ipBo, which is mapped from the 
ground state "0ideaio of the trapped ideal Fermi g as, a 
Slater determinant of the lowest N single-particle eigen- 
functions <j> n of the harmonic oscillator (HO): 

I (JV-1,JV) 

V'idealO^l, ■ ■ ' ,%n) = ~ 1= det <j) n {Xj) . (5) 
VN\ (n.j)=(0,l) 



The HO orbitals are 

4>n{x) 



ir^xlilVZ^. 



e-Q 2 / 2 H n (Q) 



(6) 



with H n (Q) the Hermite polynomials and Q = x/x osc . 
The ground state is a van der Monde determinant [2(| 

(N-l,N) (JV-l.JV) 

det H n ( Xj ) = 2 N{ - N -^' 2 det (x,) n 
(n,j) = (0,l) (f»,j-)=(o,i) 

= i n(n - i)/2 n (^-^),(7) 

yielding an exact analytical expression for the A^-boson 
ground state: 



i>Bo(xi,ar, • • • ; x N: a N ) = C_ 



N 

N n« 

.2=1 

x J~[ a(xj,<Tj;xt, at)(xj - xt) 

l<j<£<N 



(8) 



with normalization constant 
C N = 2 N ^ N - 1 ^ i ( — ) 

\ Xosc J 



N-l 



-1/2 



(9) 



This is an exact ground state in the TG limit A = +00. 
It is one member of the S^-fold degenerate ground man- 
ifold, and other members are generated by changing the 
signs of various terms in a, as in the case N — 2. If the 
product defining M were over nonover lapping pairs (1,2), 
(3, 4), • • • (N ~ 1, N) then this would give (S 2 )^ 2 = 3 N 
linearly independent states, covering the ground mani- 
fold, but pair overlap in the product is a complication; 
however, it seems reasonable to conclude that the ground 
manifold is covered. 

Suppose that now A is changed instantaneously to —00. 
Then the state © is still an exact energy eigenstate, since 
H commutes with M except at like-spin collision points 
Xj — xi, where the wave function vanishes with a cusp. 
In fact, this is an exact energy eigenstate in the limit 
|A| — > 00 even in the dissipative case where A is complex, 
since the wave function vanishes. For A <C — 1 such a 
state is highly excited, the much lower ground state b eing 
a hybrid of a totally collapsed McGuire cluster state 12\ 
for like spins and an ideal Fermi gas ground state for 
unlike spins, as in 16 1. 

Suppose that A is switched from large positive to large 
negative values rapidly enough that the sudden approx- 
imation is valid. Then the initial wave function will be 
nearly equal to the state of Eqs. ([3]) and flSJ), and the 
wave function after the switch will be a superposition of 
all eigenstates of the system with the given negative A 
value, Mt) = E a (^\ipBo)ip\ a e- iE ^ h . The domi- 
nant term in ip\ will be an STG state which reduces to 
([S]) as A — s- —00. The obvious generalization of the N=2 
STG wave function to arbitrary N differs from Eqs. ^ 
and ([8j) through replacement of \xj — xg\ by the N = 2 
solution D v for equal spins, while retaining the mapping 
from the ideal Fermi gas form for unequal spins: 



where 



n 

,aj;xe, at) 

i<i<^<w 




2 



(10) 



P v (xj,aj;x e ,ai) = ^ $<j lS 5<7 tS D u (\qje\)e q2 > e/ ' i 

s=-l 

+ ( S >? j s 1 5a e s 2 - 5 ajS2 8 aiSl ){x,j - X t ) (11) 

-1<S1<S 2 <1 

and qji = (x,j — Xi)/x osc . It satisfies the contact con- 
ditions exactly. One expects the existence of a highly 
excited gas-like STG state with nodes only at a near- 
est neighbor separation \xji\ = x no d e when the spins 
are equal, as in the spinless Bose case 14[, where a; no dc 
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increases with decreasing |A|, is very close to a\o for 
|A| 3> 1, and goes to zero along with auj in the TG limit 
|A| — > oo. For all N > 2 the approximate wave func- 
tions (p~0|) have exactly these properties, vanishing only 
at \xji\ — a: no de when the spins are equal, and becoming 
exact both at the collision points Xji = and when all 
\xjg\ — > oo. Hence we expect the unknown exact STG ex- 
cited state for finite negative A to be well approximated 
by the state of Eq. (flQ|) . In the exterior region where all 
\xjg\ > :r no dc, the wave function of this state is identical 
with that of a gas of hard spheres of diameter x no d e apart 
from normalization, and the energies of these two states 
are exactly equal. The proof is the same as those given 
previously for the STG states of spinless bosons [ljj and 
spinor fermions fljj ]. 

There also exists a quite different STG state when 
= iVj, = N/2. In this case it was shown recently 
[la ] , for a system on a ring with no circumferential trap- 
ping potential, that for — oo < A <C — 1 there exists a 
gas-like energy eigenstate in which all the fermions are 
tightly bound into (^1) dimers which behave as strongly 
attracting bosons, forming an exact analog of the un- 
trapped solution for spinless bosons 22J. It is a gas- 



[4] 

[5] 
[6] 



like STG state lying much lower than the fermionic STG 
state of Eq. (fTQ|) due to the negative binding energy of 
the dimers, but much higher than a McGuire-ideal Fermi 
hybrid state similar to that of [lq . One expects that a 
similar STG state exists for the present case of harmonic 
trapping. In fact, for N — 2 it has already been given as 
the Nf = N± = 1 case of the previously-discussed exact 
N = 2 solution, and for larger even N there will be a 
Bose-like STG state similar to that of 22[. However, it 
could not be created by a sudden switch from a strongly 
repulsive to strongly attractive interaction as in (l3| . In 
fact, the sudden approximation probability of finding the 
system in such a state after a switch A^>1^— A^;— 1 
is exactly zero in the limit |A| — > oo, since the state (fT0|) 
is an exact energy eigenstate for both A = +oo and 
A = — oo. Then by continuity this probability will be 
•C 1 after a sudden switch A ^> 1 — > —A <C — 1. 
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